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ABSTRACT 

In this work, we raise the hypothesis that the density fluctuations field 
which originates the growth of large scale structures is a combination of two 
or more distributions, instead of assuming the widely accepted idea that the 
observed distribution of matter stems from a single Gaussian field produced in 
the very early universe. By applying the statistical analysis of finite mixture 
distributions to a specific combination of Gaussian plus non-Gaussian random 
fields, we studied the case where just a small departure from Gaussianity is 
allowed. Our results suggest that even a very small level of non- Gaussianity 
may introduce significant changes in the cluster abundance evolution rate. 



Subject headings: Cosmology: large-scale structure of the universe - cosmology: 
theory — galaxies: clusters 



1. Introduction 

Generally, the problem of structure formation is associated to the gravitational growth 
of small density fluctuations generated by physical processes in the very early universe. 
Also, these fluctuations are supposed to build a Gaussian random field (GRF), where 
the Fourier components 5^ have independent, random and uniformly distributed phases. 
Such a condition means that phases are non-correlated in space and assures the statistical 
properties of the GRF are completely specified by the two-point correlation function or, 
equivalently, by the power spectrum P(k) = \5k\ 2 , which contains information on the 
density fluctuation amplitude of each scale k. This makes the choice of a GRF the simplest 
initial condition for structure formation studies from the mathematical point of view. 
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At the same time, the GRF simplicity is vindicated by a great number of inflationary 
models that predict a nearly scale-invariant spectrum of Gaussian density perturbations 
from quantum-mechanical fluctuations in the field that drives inflation (Guth & Pi 1982). 
Likewise, the central limit theorem guarantees a GRF if a wide range of random physical 
processes acts on the distribution of matter in the early universe. 

However, a number of mechanisms can generate non-Gaussian density fluctuations. 
For instance, they arise in some inflation models with multiple scalar fields (e.g. Salopek, 
Bond & Bardeen 1989); or after phase transitions when different types of topological defects 
can be formed (Kibble 1976); still, by any discrete, random distributed seed masses like 
primordial black holes and soliton objects (Sherrer & Bertschinger 1991); as well as in 
astrophysical processes during the non-linear regime where early generations of massive 
stars produce shocks which sweep material on to giant blast waves triggering formation 
of large-scale structure (Ostriker & Cowie 1981). Thus, in order to better understand 
the process of structure formation, it is necessary to investigate the possibility of the 
non-Gaussian statistics contribution to the density fluctuation field. 

Due to the difficulty to work with generic statistical models, the usual approach is to 
examine specific classes of non-Gaussian distributions. Examples of these efforts are the 
studies carried out by Weinberg & Cole (1992) that studied non-Gaussian initial conditions 
generated by a range of specific local transformations of an underlying Gaussian field; 
Moscardini et al. (1991) investigated whether non-Gaussian initial conditions can help to 
reconcile the CDM models with observations; and Kayama, Soda & Taruya (1999), who 
used data on the abundance of clusters at three different redshifts to establish constraints 
on structure formation models based on chi-squared non-Gaussian fluctuations generated 
during inflation. 

In this work, we propose a new approach to this problem, exploring the hypothesis that 
initial conditions for structure formation do not build a single GRF, but a combination of 
different fields, produced by different physical mechanisms, whose resultant effect presents 
an arbitrarily small departure from the strict Gaussianity. The paper is organized as 
follows: in Section 2, we introduce the statistical analysis of finite mixture distributions and 
present a two-component mixture model; in Sections 3, we apply the model to the cluster 
abundance evolution; in Section 4 we summarize and discuss our results. 
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2. Mixture Distributions Models: The Positive Skewness Case 

Suppose the density fluctuations field, given by the density contrast 5 = (p(r) —~p)/~p, 
is a random variable which takes values in a sample space dt, and that its distribution can 
be represented by a probability density function of the form 

p(5) = (*) + ... + «*/*(<*) (5 eft) (2.1) 

where 

aj > 0, j = l,...,k; ai + ... + a k — 1 

and 

fj(S)>0, I f,{5)d5 = 1, j = 1, 

When this happens, we say that 5 has a finite mixture distribution defined by (2.1), where 
the components of the mixture are fi(S), fk(8) and the mixing weights are aii, otk (e.g. 
Titterrington, Smith & Makov 1985). Note that we are not using here the central limit 
theorem. Mathematically, this will be valid only when k — > oo and the weights have similar 
values, so that one process has no more importance than the others. We are not making 
these hypotheses here and, consequently, the summation of processes will not necessarily 
converge to a Gaussian. 

Statistical evidence for a small level of non-Gaussianity in the anisotropy of the 
cosmic background radiation temperature has been found in the COBE 4 year maps (e.g. 
Ferreira, Magueijo & Gosrki 1998; Pando, Valls-Gabaud & Fang 1998; Magueijo 1999). 
Non-Gaussian statistics is also expected in the framework of biased models of galaxy 
formation (Bardeen et al. 1986). In this case, analytical arguments show that non-Gaussian 
behaviour corresponds to a threshold effect superimposed on the Gaussian background 
(Politzer & Wise 1984; Jensen & Szalay 1986). In the same way, hybrid models show that it 
is possible for structure to be seeded by a weighted combination of adiabatic perturbations 
produced during inflation and active isocurvature pertubations produced by topological 
defects generated at the end of the inflationary epoch (e.g. Battye & Weller 1998). Thus, a 
very compelling way to simplify our model is to apply (2.1) to the combination of only two 
fields: a GRF plus a second field, where the latter will represent a small departure from the 
strict Gaussianity. This can be posed as 

p(5) = a/i(*) + (l-a)/ 2 0$) (2.2) 

The first field will be always the Gaussian component and a possible effect of the second 
component is to modify the GRF to have positive and/or negative tails. The parameter 
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a in (2.2) allow us to modulate the relative importance between the two components. It 
represents an arbitrarily small departure from the strict Gaussianity and can be due to 
some primordial mechanism acting on the energy distribution. Such a two-component 
random field can be generated by taking 5% = P(k)u 2 , where v is a random number with 
distribution given by (2.2). Then we have 



(6 2 (r)) 



V 



P(k) 



[ah{v) + (1 - a)f 2 (v)]v 2 dis 



d 3 k 



(2.3) 



where V is the volume of an arbitrarily large region of the universe. The quantity in the 
brackets will be defined as the mixture term 



r**= f[ah{v) + {l-a)f 2 {v)\v 2 dv 



(2.4) 



so that P{k) m%x = P(k)T mtx , for the case where a is not scale-dependent. In the same way, 
the rms mass overdensity within a certain scale R will be a 2 (R) mix = a 2 (R)T mix . 

As an ilustration, in this work we explore the case of a positive skewness model, where 
the second field adds to the Gaussian component a positive tail representing a number of 
rare peaks in the density fluctuation field. A simple way to obtain this effect is to take 
the well known lognormal distribution as the second component. Besides its mathematical 
simplicity, this distribution seems to play an important role over the non-linear regime for 
a wide range of physical scales (e.g. Coles & Jones 1991, Plionis & Vardarini 1995, Bi & 
Davidsen 1997). Accordingly, our mixture becomes 



hiy) 
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e 2 



(2.5) 



(for the case of mean zero). Introducing (2.5) in (2.4), we find 
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(2.6) 



Resolving this integral we have 



rpmiX 



a+—(l-a) 



(2.7) 



Hence, if a 1, then P(k) m,tx P{k) and a 2 (R) mix o- 2 (R), which means that a 
sufficiently small contribution of the second field leaves the amplitude and shape of the 
power spectrum and the mass fluctuation practically unchanged. 
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3. Cluster Abundance Evolution 



The correct framework to describe the evolution of non-linear objects in the context of 
this model requires a generalization of the Press & Schechter formalism (Press & Schechter 
1974) in order to take into account the second field. Assuming that only regions with v > v c 
will form gravitationally bound objects with mass larger than M by the time t, the fraction 
of these objects can be calculated through 

roo 

F(M) = / p(u)du, (3.1) 



where v = S/a R . This quantity is transformed into the comoving number density of objects 
with mass between M and M + dM by taking dF/dM and dividing it by (M/p b ). Thus, 



n (M)dM = 2(^-)^- 



roo 

/ p{y)dv 

JVr. 



dM 



(3.2) 



where p^ is the background density and the number 2 comes from the correction factor 
[/o°°P( z/ )^ z/ ]~ 1 = 2, which takes into account all the mass of the universe. If p{v) is given by 
(2.2), then (3.2) can be written as 



n(M)dM = 2(^)^- 
Now, introducing (2.5) in (3.3) we have 



/ [a/iW + (1 - a)f 2 {y)]dv 
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(3.3) 



(3.4) 



Following Sasaki (1994), we rewrite (3.4) to give the density of objects with mass in the 
range dM about M which virialize at the redshift z and survive until the present epoch 
without merging with other systems. It becomes 



n(M,z) = F(Q) 



M 



(n + 3) 
3 



vf (ih) M(M? z) + (1 " a)B(M > z)] (3 - 5) 



where 



F(Q) 



-n 



(i + l") 
(i + |n + fn^ 
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(n+3) 
6 



2 



B(M, z) = exp —-In 



( 



M 



) 



and M*{z) = M*(l + 2 )- 6 /("+ 3 ) 



Eq.(3.5) allows us to compare the cluster abundance evolution with observational 
data. Clusters, as the most massive collapsed structures, correspond to rare peaks in 
the primordial density field and so their abundance is sensitive to the occurence of 
non-Gaussianity in the density fluctuation distribution. Also, cluster evolution provides a 
constraint on the amplitude of the mass fluctuation at 8 h' 1 Mpc scale, a$, and on the 
cosmological density parameter, Q m , through the relation cr s Q^ ~ 0.5 (e.g. Henry & 
Arnaud 1991; Pen 1998). In a recent work, Bahcall (1999) shows that several independent 
methods based on clusters data indicate a low mass densitiy in the universe, Q m ~ 0.2 and, 
in consequence, cr 8 ~ 1.2, breaking the degeneracy between these parameters. 

Here, we compare the behaviour of the cluster abundance evolution given by (3.5) with 
data compiled by Bahcall & Fan (1998). As an example, we plot in Figure la some fits to 
the observational data for two different values of Q m (0.2 and 1.0). Note that our model is 
very sensitive to the parameter a. Even for (1 — a) ~ 10~ 3 — 10~ 4 (i.e., almost Gaussian 
initial conditions), the curves diverge significantly from the strict Gaussian cases. This 
means that even very small deviations from Gaussianity may introduce a significant change 
in the cluster abundance. Actually, the presence of the second field tends to slow down the 
cluster abundance evolution at high redshifts. In the case of Q m = 1.0 this effect is dramatic 
for z > 0.3, while in the case of Vl m = 0.2 the difference is less pronounced and it is clearer 
for z ^ 0.6. Indeed, by plotting the 68% confidence limits around the curve Vl m = 0.2, we 
see that Gaussian and non-Gaussian models are not clearly distinguishable for z < 1 (see 
Figure lb). This is associated to the small number of observational points and, possibly, to 
the simplicity of our model. However, even considering these caveats, our results seem to 
indicate that small deviations from the strict Gaussianity may play an important role in 
the cluster abundance evolution. 



We presented the first results of a study concerning small deviations from Gaussianity 
in the primordial density field. Using very simple arguments, we developed a model based 
on the combination of two random fields in order to take into account the non-Gaussianity 
effects. This model is physically motivated in the context of hybrid models, as well as in 
the framework of biased scenarios for structure formation. The weighted combined field 



4. Summary and Discussion 
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involves a parameter a which modulates the relative importance of its components. For 
a « 1, we preserve the amplitude and shape of -P(fc) and cr(i?) almost the same as in the 
Gaussian case. At the same time, our results suggest that even very small values of (1 — a) 
can introduce a significant change in the cluster abundance evolution. This effect seems to 
be stronger in high density universes (at z < 1) than in low density universes where the 
effect probably turns more important at higher redshifts. 

The model has some drawbacks. Firstly, it depends on the choice and amplitude of 
the second component of the combined field. Our choice of the lognormal function had 
a mathematical criterion of simplicity. A detailed investigation of the use of different 
distribution functions as the second component will be the subject of future works. However, 
the reasonable agreement between the model and the data gives some support to our 
arbitrary choice. Other possible limitation of this work comes from the use of the analytical 
approximation to the density of non-linear objects following Sasaki (1994). A more accurate 
description of the cluster abundance evolution requires the utilization of numerical methods. 
But Blain & Longair (1993), also working in the Press & Schecter framework, found results 
numerically similar to Sasaki's, so it seems that using this analytical approximation does 
not introduce any systematical error. Finally, we should keep in mind that our results are 
preliminary and both theoretical and observational efforts are necessary in order to confirm 
or disproof the hypothesis that the primordial density field can be described as a slightly 
non-Gaussian distribution. 
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Fig. la - Galaxy cluster abundance in the two-component model for Q m = 0.2 with 
a = 0.9990 (solid) and a = 1 (dashed) and for Vl m = 1.0 with a = 0.9999 (solid) and a = 1 
(dashed). The curves, normalized at z — 0, correspond to n — —1.0, M* = 10 14 h~ l Mq 
and M > 8 x 10 14 h" 1 Mq. The observational points were taken from Bahcall & Fan 
(1998). 
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Fig. lb - Galaxy cluster abundance in the two-component model for Vl m = 0.2 with 
a = 0.9990 (solid) and a = 1 (dashed). The dotted lines are the 68% confidence limits 
around the non-Gaussian fit. The observational points were taken from Bahcall & Fan 
(1998). 



